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Abstract. We present an interpolative method for describing colierent transport through 
an interacting quantum dot. The idea of the method is to construct an approximate 
electron self-energy which becomes exact both in the limits of weak and strong coupling 
to the leads. The validity of the approximation is first checked for the case of a single 
(spin-degenerate) dot level. A generalization to the multilevel case is then discussed. 
We present results both for the density of states and the temperature dependent linear 
conductance showing the transition from the Kondo to the Coulomb blockade regime. 



1 Introduction 

The Kondo effect constitutes a prototypical correlation effect in condensed mat- 
ter physics. Although originally studied in connection to magnetic impurities 
in metals, there is now a renewed interest in this many-body problem fostered 
by the recent observation of Kondo effect in semiconducting quantum dots [1,2]. 
Quantum dots provide an almost ideal laboratory where the relevant parameters 
can be controlled, which allow to test the predictions of theoretical models. 

From the theoretical side, Kondo physics in quantum dots has been mainly 
analyzed in the light of the so called single level Anderson model. There were 
predictions for the Kondo effect in quantum dots based on this model since the 
early 90's [3,4]. The theory predicts an enhancement of the linear conductance 
due to Kondo effect at very low temperatures, which is in qualitative agreement 
with recent experiments. 

However, in most realistic situations, the single level Anderson model consti- 
tutes a crude approximation for a quantum dot. Actual semiconducting quan- 
tum dots contain a large number (~ 100) of electrons and the single-particle 
level separation between dot levels may be not so large compared to the level 
broadening, which restricts the validity of the single-level approximation. The 
actual situation would be more appropriately described by a multilevel model, 
including several instead of a single dot level. Unfortunately, there are no simple 
theoretical approaches to extract the electronic and transport properties from 
such a microscopic model. 

In this paper we present results on the Kondo effect in quantum dots based 
on the interpolative method. The basic idea of this method is to construct an 
interpolative electron self-energy which becomes exact both in the limits of weak 
and strong coupling to the leads. These ideas were first introduced in Ref. [5] 
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in connection to the single-level Anderson model and have, since then, been 
adapted by several authors to different problems involving strongly correlated 
electrons. In this way, the method has been used to study the Hubbard model [6], 
the non-equilibrium Anderson model [4] , the metal- insulator transition in infinite 
dimensions [7] , to incorporate correlation effects into band-structure calculations 
[8], the ac-Kondo effect in quantum dots [9] and finally extended by the present 
authors to analyze the multilevel Anderson model [10]. 

The paper will be organized as follows: In section 2 we present the interpola- 
tive method. We first discuss the single level case, showing the accuracy of the 
method with the help of a simple exactly solvable model. We then consider the 
multilevel situation. In section 3 we present results which illustrate the behavior 
of the conductance with temperature in a multilevel situation. 

2 The interpolative method 

For describing a multilevel quantum dot (QD) we consider a model Hamiltonian 
^ = -^dot + ^leads + ^^ "^^^^^ ^dot = Em em4i'^m + UJ2i>rnf^rnni Cor- 
responds to the uncoupled QD (n™ = dl^dm); ^leads ^ ^keL,R^kCkCk to the 
uncoupled leads, and Ht = ^m keL R^'m,kdtiCk + h.c. describes the coupling 
between the dot and the leads. The labels m and / in H denote the different dot 
levels including spin quantum numbers. The number of dot levels will be denoted 
by M (i.e. 1 < m,l < M). We adopt the usual simplifying assumption of having 
the same electron-electron interaction U between any pair of dot states. 

The main objective of our method is to determine the dot retarded Green 
functions Gmir) = —i9{T) < [(i„i(T), (ij„(0)]+ > from which the different level 
charges and the dot linear conductance can be obtained. In the frequency rep- 
resentation we can write Gm as: 

Gm\^) = Tfp ^r; -, ^ — -, ^ — 7 T, (J-j 

^ - Cm - ^m[UJ) - 1 ,n.L[UJ) - 1 mM^) 

where e^^ — Cm + U Ei=:tm "' i'' ^^^ Hartree-Fock level (we adopt the notation 
ni for the mean charge on level /) and r,m,L, Tm^R are tunneling rates coupling 
the dot to the leads, given by rm,L(fl)(u;) = Y.keL(R) ^li^k/i^ - ^k + J0+). We 
shall neglect indirect coupling between dot levels through the leads (non-diagonal 
elements r^m^LiR)) and adopt the usual approximation of considering r„iHJi-j 
as a pure imaginary constant independent of the energy. 

The self-energy Z'm(a;) takes into account electron correlation effects beyond 
the Hartree approximation. The idea of the present approximation is to deter- 
mine an interpolative self-energy which yields the correct exact results both in 
the r/U -^ limit (atomic limit) and in the opposite U/F -^ limit. 

2.1 The single-level case 

Let us first discuss how to proceed for the simple single-level case. In this case 
TO = 1, 2, the two indexes corresponding to up and down spin orientations. These 
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will be denoted by a and a . In the atomic limit Ga can be obtained using the 
equation of motion technique [11] as 

This expression can be formally written in the usual Fermi liquid form, i.e. 

Ga ((jj) = [oj — e — Una — Sa (lu)]^^ by introducing the "atomic" self-energy 



uj-e-U{l-na)+iO+ 

In the opposite limit, U/F -^ 0, the electron self-energy can be calculated by 
second-order perturbation theory in U, which yields 

4^)(.) ^ u^ r de, r de. r de, ^-(-^)/^^(-^)^^(-3) ^ 

J-oo J-oo J-oo UJ + e2-ei-€'i+lQ+ 

[/(ei)/(e3) (1 - /(e2)) + (1 - fie,)) (1 - f{e,)) f{e,)] , (3) 

where f{uj) is the Fermi distribution function and Pa{^) = r /'k{{uj — laY + r^) 
is the local density of states for an effective level eg, which will be determined 
in order to fulfill exact Fermi liquid properties at zero temperature. 
It is important to stress the following simple property of S^"^': 



r^o '^ ' ' uj -ea + iQ+ 

Thus, when extrapolated to the atomic limit S^^> has the same functional 
form as S^°-*\ This property suggests that one can smoothly interpolate between 
the two limits. The ansazt proposed in Ref. [5] for the interpolative self-energy 
is: 

1 - a Sa (uJ) 

where a — (e — ia — U{1 — ng-))/(f7^ng-(l — n^))- This ansazt has the desired 
property S -^ Z'^^) when f7 ^ and Z" ^ U^"*^ when T ^ 0. 

The final step is to impose the proper self-consistent condition for determin- 
ing the effective level e. At zero temperature, from the Luttinger-Ward relations 
[13] one can derive the Friedel sum rule for the Anderson model [14] 

n,, = --Im In G;(^f) 

TT 

which imposes an exact relation between the dot-level charge and the phase 
shift at the Fermi energy. The effective level can thus be determined in order 
to fulfill the Friedel sum rule. This condition is, however, not valid at finite 
temperature. In Ref. [4] we show that the condition Ua — fia, i.e. imposing the 
same charge in the effective system as in the interacting system, is approximately 
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equivalent to the Fricdcl sum rule at zero temperature but can be also used at 
finite temperature. 

In order to check the accuracy of the interpolative method we have considered 
a simple two-sites problem that can be diagonalized exactly. One of the sites 
would describe the metallic leads and the other site corresponds to the dot. In 
order to analyze the more general situation we impose a finite splitting A = 
ta — eg. between the two spin orientations on the dot. Within this toy model the 
second order self-energy can be evaluated analytically. 

In figure 1 we show the charge on the two dot levels as a function of gate 
voltage (the gate voltage is the distance between the lower dot level and the 
leads level). As can be observed, in the exact solution there is a blocking of the 
upper level charge until the gate voltage becomes larger than A + U. The exact 
behavior is accurately reproduced by the interpolative method. It is instructive 
to consider another simple approximation widely used in the literature, which 
consist in just broadening the poles in the atomic Green function (2) by the 
non-interacting tunneling rates. This approximation corresponds to the so-called 
Hubbard I approximation [12]. As can be observed in the lower panel of Fig. 1, 
this approximation fails to give the blocking of the upper level found in the exact 
solution. 
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Fig. 1. Level charges as a function of gate voltage for the two sites model with A — 0.25 
and i = 0.1 (in units of the charging energy U). Left pannel corresponds to the inter- 
polative approach and the right pannel to the Hubbard I approach. The exact solution 
is shown as a dashed line. 



2.2 Multilevel case 



The multilevel version of the interpolative method is somewhat more complex 
[10]. In the first place, the atomic limit Green functions do not contain just two 
poles but several poles corresponding to the various different charge states of the 
dot. The corresponding expression can be obtained using the equation of motion 
technique and is given by 
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OJ — €jn ^ ^0+ f-^ UJ — €jn - U ^ 20+ 



LU-era- {M -1)U + iO+' 

The evaluation of this expression requires the knowledge of up to M — 1-body 
correlations functions < hin2 >, < nin2n3 >, ..., etc. However, for sufficiently 
large U fluctuactions in the dot charge by more of one electron with respect to 
the mean charge TV are strongly inhibited. One can thus approximate Eq. (5) as 
follows 



Gtf)(^) 



Am Am 



UJ-Ern- U{N -l)+iO+ UJ- £„ -UN + iO+ 

Am 

^^+1 (6) 



U{N + l)+iO+' 



where N = Int[M]. In order to yield the exact first three momenta of the exact 
spectral density the weight factors A^ should satisfy the following sum rules 



Am I Am I Am _ i 

{N - 1)^A^_, + N^A^ + {N + lfA^+, = I] "i+ < "^ >m, (7) 

l^m 

where < nn >m= ^n-ik)^m < f^ifik >■ For the special case N — {N = M — 1) 
one has ^™_i = (^jv+i — 0) and only the first two Eqs. in (7) have to be 
considered. 

This approximated expression for Gm {^) is now fully determined by the 
average charges ni and the two-body correlation functions < hifik >■ As in 
the single-level case one can define an atomic self-energy, Sm ^ uj — e^f — 



Gt'\Lo) \ 



which can be written as the ratio of two polynomials in lj 



^{at) ^ amU'^{uJ-em+lO+)+hmU^ ,^. 

[uj — e„i + iO+)2 + CmU{uj — €„L + «0+) + d„iU^ ' 

where a™ = (A/" - «,„) [1 - (A/" - n^)] + < hh >,„; c^ == TV - n^ - 3TV; d.^ =< 
nh>m +3TV2- 1- (3iV- l)(TV-n„) and 6„ ^ N^{1- N) - (Af - nm)dm- 

On the other hand, in the U/ F — > limit the self-energy is accurately given 
by second order perturbation theory as in the single level case. The second order 
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(2) 

self-energy Sm now takes into account the interaction of an electron on the dot 
level m with electron-hole pairs on each one of the other channels. 

For the interpolation one notices that both Z"'^^ and Z'^"*) have the same 
functional form when extrapolated to the corresponding opposite limit. The 
natural generalization of the ansazt in the single level case now has the form of 
a continued fraction 

2^m[l^) — JY] 7 nm(^) 



with coefficients a™ = U'^am/Am, Pm = (e™ - ?m + {bm/am - Cm)U)/A. 



■m-i 



7m = {{cyyi - b„i/am)bm/am - dm)U / A^ and dm = (em - e„i + Ubm/arn)/Am, 
where A^ ^ U^ Y^i^m "-'(1 - "O- 

As in the single-level case the final step is to determine the effective levels 
self-consistently. In the multilevel case one has, in addition to self-consistently 
determine the two-body correlations < iin >m that appear in the atomic self- 
energy. This can be done by means of the relation 

1 1"°° 
V" < fiifim >= —I .f (t^)Im [{uj -era- r,n) G„(a))] dw, (9) 

connecting the two-body correlations and the Green functions that can be de- 
rived from the equation of motion for Gm(w). This step turns out to be essential 
in order to obtain the correct values of the charges in the large U limit. 

Finally, the temperature dependent dot linear conductance can be obtained 
using the expression [15] 

G=-r2_^ 7T7^ HTTT^ h / [tTI ImG„ w)dw 



3 Results 

The multilevel formalism allows to study the importance of the multilevel struc- 
ture in the QD transport properties. For this purpose we have analyzed the 
M = A case which corresponds to two consecutive dot levels plus spin degener- 
acy. We have studied this case as a function of the level separation A. 

Figure 2 shows the dot conductance as a function of Fermi energy and tem- 
perature for the cases Z\ = 0, 0.1 and 0.5 (in units of U). 

This figure illustrates the transition from a two-fold degenerate situation 
(Z\ = 0), where the conductance reaches a maximum value 4e^//i for the half- 
filled case at zero temperature, to the case of well separated dot levels, where 
the maximum conductance 2e^//i is reached for the quarter and three quarter 
filling case. The increase of conductance with decreasing temperature is due to 
the Kondo effect. While in the case of well separated levels one observes only 
the Kondo effect due to the spin-degeneracy of the individual levels, when A is 
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Fig. 2. Dot conductance as a function of Fermi energy for increasing temperatures 
with Fl ~ Fr — 0.075 and Zi = (a), 0.1 (b) and 0.5 (c). The temperature values are 
0.0005, 0.0025, 0.005 and 0.03 in units of U 



small compared to F one can observe Kondo features involving the two nearby 
dot levels. When the temperature is raised above the Kondo temperature (which 
is around 0.005 for the parameters used in this figure) one recovers the sequence 
of dot resonances at the charge degeneracy points characteristic of the Coulomb 
blockade regime. 

The Kondo effect should manifest also as a zero-bias anomaly in the dot 
non-linear conductance. This anomaly is directly related to the appearance of a 
narrow peak around the Fermi energy in the dot spectral density. In cases where 
the splitting between dot levels is of the order of F we expect to have a zero-bias 
anomaly not only between dot resonances corresponding to the same dot level 
but also in between resonances corresponding to different levels. This feature is 
illustrated in Fig. 3 where we plot the density of states around Ep for the same 
three cases of Fig. 2 with Ep = 1.5. The appearance of a zero-bias anomaly in 
between resonances corresponding to different levels is a clear manifestation of 
the multilevel structure of the QD which has already been observed in recent 
experiments on semiconducting quantum dots [16]. 

This work has been funded by the Spanish CICyT under contracts PB97-0028 
and PB97-0044. 
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Fig. 3. Density of states around the Fermi energy for increasing temperatures values 
(0.01, 0.02 and 0.03) for the same three cases in Fig. 2 and Ef ~ 1.5 
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